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"An  analytical  nodal  for  calculating  transfer  functions  relating  the  C-130 
lightning  sensor  responses  to  incident  electromagnetic  fields  was  developed. 
The  model  consists  of  two  components;  a  low  frequency  non-resonant  term  and  a 
frequency-dependant  tarn  including  airframe  resonances.  The  latter  are 
calculated  for  the  frequency  range  1-10  >Ss  and  compared  with  measured  data. 
Recommendations  for  removing  airframe  resonances  from  the  measured  data  are 
presented. 
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This  final  technical  report  was  submitted  by  the  Mechanical /Electrical  Systems 
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Brian  Ruhlman,  AFWAL/F1ESL,  as  the  Project  Engineer.  W.  P.  Geren  was  the 
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SECTION  I 
INTRODUCTION 

This  report  describes  an  analytical  model  for  the  calculation  of  total 
electric  and  magnetic  fields  induced  on  the  surface  of  an  aircraft  by  an 
incident  plane  wave  field.  Transfer  functions  are  presented  for  nine 
lightning  sensors  in  the  form  of  equations  for  the  non-resonant  component. 

For  three  sensors,  transfer  functions  are  computed  over  the  1-iO  MHZ  range  and 
compared  with  data. 

1 .  BACKGROUND 

Data  on  the  radiated  electromagnetic  fields  produced  by  lightning  strikes 
was  gathered  during  the  sumawrs  of  1981-82  in  a  joint  National  Oceanographic 
and  Atmospheric  Administration  (NOAA)  and  AFVAL/FIESL  project.  A  WC-130 
aircraft  was  equipped  with  electrically  small  antennas  (lightning  sensors)  to 
measure  surface  fields  on  the  airframe. 

As  part  of  this  program,  the  Boeing  Company  was  contracted  to  provide 
technical  assistance  In  the  calibration  of  the  lightning  sensora  by 
calculating  transfer  functions  relating  surface  fields  on  the  airframe  to  an 
incident  plane  wave  field. 

2.  OBJECTIVES 

\ 

The  overall  objective  of  the  effort  described  In  this  report  was  to 
provide  guidelines  and  quantitative  methods  for  processing  the  C-I30  lightning 
sensor  data.  These  methods  are  based  upon  the  results  of  the  analytical  model 
and  will  be  used  to  remove  the  response  of  the  airframe  from  the  sensor  data, 
resulting  in  estimates  of  the  incident  electromagnetic  fields.  Specifically 
the  objectives  were  to: 

a.  Develop  an  analytical  model  for  calculating  electromagnetic 
fields  on  the  airframe  surface. 

b.  Calculate  transfer  functions  for  all  lightning  sensors. 


Die  the  lightning  tensor  date  from  the  1981-82  flight  program 
along  with  the  corresponding  transfer  functions,  to  infer 
quantitative  methods  for  removing  airframe  effects  from  the 
data. 
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SECTION  II 
TECHNICAL  APPROACH 

Currents  induced  on  the  surface  of  a  relatively  thick  structure,  such  as 
an  aircraft  fuselage,  aay  be  decomposed  into  three  components: 

a.  Axial  currents  with  no  azimuth  variations  (resonant). 

b.  Axial  currents  with  azimuthal  variation  (non-resonant). 

c.  Azimuthal  currents  (non-re sonant). 

• 

As  described  in  reference  1,  currents  of  type  a.  may  be  obtained  from  a 
thin  wire  model  of  the  airframe.  The  corresponding  surface  current  density  is 
obtained  by  dividing  the  total  current  by  the  effective  clrcumf erence . 

Surface  currents  of  type  b.  or  c.  may  be  approximated  in  the  low  frequency 
limit  (l.e.,  wavelength»effective  radius),  by  expressions  for  scattering  from 
cylinders  (fuselage)  and  strips  (wings).  These  expressions  are  simple  and 
require  minimal  computation.  The  thin  wire  model,  however,  requires  the 
solution  of  an  integral  equation.  The  two  model  types  will  be  referred  to 
throughout  as  "non-re sonant"  and  "resonant”,  since  all  airframe  resonance 
effects  are  contained  in  the  thin  wire  model. 

One  method  for  solving  thln-wlre  problems  is  the  method  of  moments.  In  the 
method  of  moments,  the  airframe  structure  is  subdivided  into  equivalent  wire 
segments.  Coupling  between  the  segments  is  calculated,  resulting  in  an 
Impedance  matrix  which  is  then  used  to  relate  induced  currents  to  Incident 
fields.  This  technique  is  well  established  but  requires  large  computers  and 
elaborate  input  data.  We  opted  in  this  program  for  the  approach  described  in 
reference  2.  In  this  alternate  method,  the  integral  equation  is  solved  by 
obtaining  approximate  analytical  solutions.  In  the  lowest  (zeroth)  order 
approximation,  the  current  on  each  airframe  member  i.e. ,  forward  fuselage, 
right  wing,  left  wing,  aft  fuselage,  vertical  stabilizer,  right  horizontal 
stabilizer,  and  left  horizontal  stabilizer,  is  expressed  as  the  sum  of  a 
sinusoidal  term,  with  undetermined  coefficients,  and  a  forced  term, 
proportional  to  the  tangential  electric  field.  Imposing  current  and  charge 
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conservation  at  the  intersections,  along  with  the  boundary  condition  that  the 
current  goes  to  sero  at  wing  tips,  fuselage  nose,  and  stabilizer  tips, 
determines  the  sinusoid  coefficients.  The  first  order  correction  is  obtained 
by  the  Iteration  procedure  described  in  reference  3.  The  thin-wire  model 
developed  under  this  contract  Includes  the  first  order  correction  and  is  in 
good  agreement  with  results  published  by  R.  W.  P.  King.  The  results  shown  in 
Figs.  I  and  2  for  symmetric  and  antisymmetric  currents  induced  in  an  equi-arm 
cross  are  to  be  compared  with  Ref.  3,  Fig.  3. 


Symmetric  current  at  center  of  equi-ara  cross  calculated  from 
thin -wire  model  for  comparison  with  published  results  (cee  text) 


symmetric  current  at  center  of  equi-arm  cross  calculated 
thin-wlre  model  for  comparison  with  published  results  (see 


SECTION  III 

DESCRIPTION  OF  ANALYTICAL  MODEL 


As  described  previously,  Che  surface  currents  may  be  obtained  from  two 
independent  models,  one  which  gives  the  net  current  on  each  airframe  member, 
and  another  which  gives  azimuthally-varying  and  azimuthally-directed  currents. 
The  total  surface  current  at  a  given  location  is  the  vector  sum  of  the  two. 

The  surface  charge  density,  q,  is  obtained  from  the  conservation  law: 

divJ  +  jwq  -  0 

where  J  is  the  surface  current  density,  div  is  the  divergence  operator,  and  w 
is  the  radial  frequency. 

Choose  the  aircraft  coordinates  as  follows  (see  Figure  3): 

x-axis  directed  forward  along  fuselage 

y-axis  directed  along  the  left  wing 

z-axis  upward 

and  the  lightning  channel  coordinates  as  follows: 

x-axis  parallel  to  y-axis 

y-axis  anti-parallel  to  x-axis 

Z-axis  parallel  to  z-axis 

.  origin  at  location  of  ground  strike 

RKo,  theta,  phi  are  the  spherical  coordinates  of  the  aircraft  in  the 
(x,?,z),  or  lightning  channel  reference  system. 

1.  THIN  WIRE  MODEL 


The  thin  wire  model  provides  the  net  current  and  charge  on  the  airframe 
at  a  given  location.  The  corresponding  surface  current  and  charge  are 
obtained  by  dividing  by  the  effective  circumference.  For  the  fuselage,  which 
is  cylindrical,  the  effective  circumference  is  the  geometrical  circumference. 
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For  Che  the  wings,  Che  effective  circumference  is  Che  geometrical 
circumference,  multiplied  by  a  concentration  factor  which  varies  from 
typically  .5  in  mid-chord  to  4  or  5  near  an  edge. 

2.  CYLINDER  AND  PROLATE  SPHERICAL  MODELS  (FUSELAGE) 

Using  the  approach  described  in  reference  1 ,  the  surface  current  density 
induced  on  a  cylinder  with  axis  along  the  x-axis  is  given  by: 

J  -  <JA  +  2  3  Klnc)x  -  xHlnc  (1  -  2jk-r]0 

where:  H  c  *  incident  magnetic  field 

JA  ■  total  axial  current/circumference 

9  -  -sinG  y  +  cosG  z 

r  *  a(cos0y  +  sinGz) 

tf  ■  wave  number 

a  *  radius  of  cylinder 


The  normal  electric  field  is  obtained  from  the  conservation  of  charge  and 
is  given  by: 


E  .  +  2n  E. 
nA  Inc 


where:  ®inc  *  incident  electric  field 

EnA  "  qA/BPSO,  qA  -  (j/kc)dJA/dx 

c  ■  velocity  of  light 

EPSO  ■  permittivity  of  free  space 
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For  E-field  antennae  located  near  the  extremeties  of  the  fuselage,  the 
prolate  spheroid  model  is  of  use.  The  normal  electric  field  induced  on  the 
surface  of  a  prolate  spheroid  of  major  and  minor  axis  a  and  b,  respectively, 
is  given,  in  the  quasistatic  approximation,  by  (major  axis  along  x-axis): 


F«  +  Fr  Wr'*> 


where:  (x/a)2  +  (r/b)2  ■  1 


A 

r 


cos©  y  +  sinQ  z 


Fx,Fr»  functions  of  eccentricity,  e 


(1  -  (b/a)2)*5 


Note  that  the  quasistatic  prolate  spheroid  model  includes  coupling  to  the 
axial  component  of  Einc  which  ia  the  low  frequency  limit  of  the  resonant  model. 

3.  STRIP  MODEL  (WING) 

The  surface  currents  induced  on  the  upper  surface  of  the  wing  are 
obtained  from  the  thin  wire  model  and  the  infinite  strip  (see  reference  4)  and 
are  given  by  the  expression: 

• 

JA  +  Hxinc  "  V^inc 

total  axial  current /effective  circumference 
Incident  magentic  field 

t/((V/2)2  -  t2r5 

wing  width 

distance  from  mid-chord  to  sensor  in  the  positive 
x-direction 


where: 


inc 

*W 

W 
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SECTION  IV 

DESCRIPTION  OF  SENSOR  TRANSFER  FUNCTIONS 

The  nine  lightning  sensors  are  listed  in  Table  I.  As' indicated  in  the 
table,  the  sensors  measure  the  tine  derivatives  of  the  surface  electromagnetic 
fields.  According  to  the  vendor’s  data  sheets,  the  sensor  outputs  into  the 
specified  load  are  proportional  to  the  tine  derivative  of  the  field  at  the 
sensor,  and  ail  have  a  bandwidth  (3  dB  point)  in  excess  of  38  MHz.  Hence,  the 
transfer  function  relating  sensor  output  to  tine  derivative  of  surface  field 
is  a  frequency-independent  factor  given  in  the  vendor's  data  sheets.  The 
remainder  of  the  overall  transfer  function  relating. sensor  output  to  incident 
field  is  given  by  the  analytical  nodel. 

The  location  of  the  sensors  determines  the  choice  of  effective 
circumference  and,  in  some  cases,  the  importance  of  the  contribution  of  the 
thin  wire  model.  The  sensors  are  discussed  individually  in  the  following. 

A.  General  Description  of  Analytical  Models  for  Sensors 

U  Sensors  #1  through  #3,  #6  and  #7 

These  sensors  are  all  located  on  the  upper  centerline  of  the  fuselage. 

The  effective  circumference  for  all  is  just  the  fuselage  circumference.  Both 
thin-wire  and  cylinder  models  are  required. 

2.  Sensors  #4  and  #9 

Sensors  #4  and  #9  are  located  just  forward  of  mid-chord  and  have  a 
concentration  factor  of  approximately  .7,  obtained  from  calculating  the 
distribution  of  fields  about  the  wing.  Both  the  thin  wire  model  and  the  strip 
model  are  required  for  these  sensors. 

3.  Sensor  #5 

This  sensor  presents  a  different  sort  of  geometry,  with  the  effective 
circumference  dependent  on  the  two  principal  radii  of  curvature  of  the  wing 
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Sensor 


TABLE  I  -  C-130  LIGHTNING  CHARACTERIZATION  SENSORS 


Model 

Type 

Use 

Location 

Naoe 

EG&G  FPD-2(A) 

D-dot 

Electric  Field 

Forward  Upper  Fuselage  E-FUF 

EG&G  CML-S7A(A) 

B-dot 

Wing  Axis  H-Field 

Forward  Upper  Fuselage  J-NT 

EG&G  CML-S7A(R) 

B-dot 

Fuselage  Axis  H-Field 

Forward  Upper  Fuselage  J-WW 

EG&G  CML-7(R) 

• 

B-dot 

Wing  Surface  Current 

Left  Upper  Wing 

J-LUW 

EG&G  RPD-2B(R) 

D-dot 

Wing  Electric  Field 

Left  Wlngtlp 

E-LWT 

EG&G  HSD— S1A(R) 

0  -dot 

8 

Electric  Field 

Aft  Upper  Fuselage 

E-AUF 

EG&G  CML-S7A(R) 

B-dot 

Fuselage  Surface  Current  Aft  Upper  Fuselage 

J-AUF 

EG&G  FPD-2(A) 

D-dot 

Electric  Field 

Aft  Lower  Fuselage 

E-ALF 

EG&G  CML-7(R)  B-dot  Wing  Surface  Current 


Right  Upper  Wing 


J-RUW 


tip.  A  reasonable  first  approximation  is  that  the  effective  radius  of 
curvature  for  charge  density  is  the  geometric  mean  of  the  two  principal  radii 
of  curvature  (in  the  plane  of  the  wing  and  normal  to  the  plane  of  the  wing). 
Only  the  thin  wire  model  is  required  for  this  sensor. 

4.  Sensor  #8 

This  sensor  is  similar  in  its  complexity  to  the  wing  tip  sensor.  The 
sensor  is  located  at  the  .junction  of  the  horizontal  stabilizer  and  fuselage, 
so  that  the  total  charge  is  the  sum  of  the  charges  induced  on  the  two.  The 
effective  circumference  is  obtained  by  considering  how  this  net  charge 
distributes  over  the  junction  area.  Below  the  first  resonance  of  the 
horizontal  stabilizer,  it  will  act  as  a  capacitive  load  on  the  fuselage,  and 
the  charge  will  spread  out  approximately  uniformly  over  the  horizontal 
stabilizer.  Both  the  thin  wire  and  prolate  spheroid  models  are  required  for 
this  sensor. 

B.  Conventions  for  Units  and  Transfer  Functions 

In  order  to  organize  the  transfer  functions  and  data,  all  transfer 
functions  will  be  divided  into  two  types.  E-field  and  H-field.  The  E-field 
transfer  function  is  defined  as  the  normal  component  of  the  E-field  at  the 
sensor  location/incident  E-field  and  is  dimensionless.  The  H-field  transfer 
function  is  the  surface  current  density  at  the  sensor  location/incident 
H-field,  and  has  two  vector  components.  Since  surface  current  density  has  the 
same  units  as  H-field,  the  H-field  transfer  function  is  also  dimensionless. 

For  the  non-resonant  models,  the  vector  components  of  the  incident  field 
are  expressed  in  the  aircraft  coordinate  system  (see  Figure  3).  For  the 
resonant  model,  the  incident  field  is  the  amplitude  of  the  plane  wave,  which 
is  described  completely  by  the  propagation  vector  and  polarization  angle.  The 
equations  for  the  non-resonant  models  are  listed  in  Table  II. 
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TABLE  II  ~  EQUATIONS  FOR  NON- RE SONANT  TRANSFER  FUNCTIONS 


1 

2 

3 

4 

5 

6 

7 

8 

9 


E-FUF 

J-NT 

J-WW 

J-LUW 

E-LWT 

E-AUF 

J-AUF 

E-ALF 

J-RUW 


En  ‘2Ezinc 


J  -  2H  . 
x  yinc 


Jy  -  (1  -  2jV)Hxlnc 


xinc 


KWH 


zinc 


negligible 


E  -  2E  , 
n  zinc 


J  -  -2H  . 
z  yinc 


E  -  -P  E  .  -PE. 
n  x  xinc  z  zinc 


Jy  *  Hxinc  'Sf^zinc 


Note:  Px,  Pz,  Ky  are  described  in  Chapter  3  under  the  prolate  spheroid  and 
strip  models.  Estimated  values  are: 


Px  “  5 

P  -  2 

z 

K,  ,  -  .4 
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SECTION  V 

RESULTS  OP  ANALYTICAL  CALIBRATION  AND  COMPARISON  WITH 
DATA  IN  THE  1-10  MHZ  RANGE 

The  equations  for  the  non-resonant  transfer  functions  of  the  analytical 
model  are  shown  in  Table  II  and  are  self-explanatory.  The  thin-wire,  or 
resonant  model  is  somewhat  more  complicated.  The  thin  wire  model  has  been 
incorporated  in  a  computer  program  which  runs  on  the  PDP  11/70.  Results  for 
three  sensor  transfer  functions,  using  the  effective  circumferences  of  Table 
III,  were  obtained  for  an  aircraft  elevation  angle  45°,  relative  to  the 
lightning  coordinate  system  (see  Figure  3),  and  a  vertically  polarized  plane 
wave  (i.e.  horizontal  H-field)  and  the  following  azimuths: 

1.  Incidence  along  the  nose. 

2.  Incidence  at  45°  to  fuselage,  between  nose  and  left  wing. 

3.  Incidence  along  the  left  wing. 

4.  Incidence  along  the  tail. 

The  thin  wire  model  may  be  used  to  obtain  transfer  functions  for  the 
remaining  sensors,  as  data  becomes  available.  The  JDOT-WW  sensor  was  omitted 
because  its  resonant  response  is  negligible,  as  it  couples  primarily  to  the 
non-resonant  circumferential  fuselage  currents. 

The  complete  set  of  C-130  resonances  is  quite  complex  (see  Ref.  5,  page 
163).  Excluding  the  effects  of  the  HP  wire,  there  are  more  than  10  airframe 
resonances  in  the  3-10  MHZ  range.  There  are  no  resonances  below  3  MHZ.  The 
thin  wire  model  was  simplified  by  representing  the  aft  fuselage  by  a  single 
cylinder,  with  an  Increased  length  to  partially  account  for  the  empennage. 

The  justification  for  this  is  that  the  major  resonances  are  determined  by  the 
wing  and  gross  fuselage  structure.  I.e.,  the  addition  of  empennage  details 
will  only  add  some  fine  structure  and  result  in  slight  shifts  of  the  major 
resonances. 

The  three  sensors  which  exhibited  appreciable  airframe  resonance  effects 
and  for  which  data  were  available  are  JNT-DOT  (or  B-DOT  NT),  JLUW-DOT  (or 
J-DOT  LUW) ,  and  EALF-DOT  (or  D-DOT  ALP).  For  each  of  the  three,  transfer 


TABLE  III  -  EFFECTIVE  CIRCUMFERENCES  FOR  RESONANT  MODEL 


E-FUF  Cp 

J-NT 

J-WW 

J-LUW  2W/.7 

E-LWT  (CiC2)*5 


E-AUF  Cp 

J-AUF  Cp 

E-ALF+  C_ 

F 

J-RUW  2W/.7 


CF  -  Fuselage 

Circumference 


W  -  Wing 
Width 

CpCj  ■  Circumferences 

corresponding  to  the 
two  principal  radii 
of  curvature  of  the 
wing  tip 


+  Complicated  by  junction  of  fuselage,  horizontal  and  vertical 


stabilizers.  CF  is  an  estimate 


functions  were  calculated  and  nultiplied  by  an  assuaed  lightning  spectrum  of 
2 

1/F  ,  to  produce  an  estimate  of  the  spectral  response  In  the  1-10  MHZ  range, 
and  to  show  the  relative  magnitude  of  the  various  resonant  peaks.  The  results 
are  discussed  In  the  following  text. 

o  JNT-DOT 

The  thin  wire  model  predicts  a  pair  of  resonances  at  3.5-4  MHZ  and 

4.2- 4. 8  MHZ  with  a  null  between  (see  Figures  4A-4D).  The  null  Is  due  to  the 

aft  fuselage-wing  resonance  and  is  described  by  R.H.P.  Ring  in  Reference  3 

(see  Figure  4  and  accompanying  discussion  on  page  515  of  Reference  3).  There 

is  a  higher  order  resonance  near  10  MHz.  The  pair  of  resonances,  and  the  null 

location  are  quite  sensitive  to  change  in  azimuth.  Comparison  with  data  (see 

Figures  7A-7C)  shows  general  agreement  with  the  lower  resonances  but  no 

pronounced  higher  order  resonance  above  6  MHZ.  Such  resonances  should  occur 

(see  Reference  1,  Figure  18)  for  an  aircraft  of  the  C-130  geometry.  Their 

absence  indicates  either  a  deficiency  in  the  lightning  spectrum  at  these 

2 

frequencies  (i.e.,  a  faster  than  1/F  roll-off)  or  inadvertent  filtering  in 
the  data  acquisition  and  processing. 

o  JLUW-DOT 

The  thln-wire  model  results  (Figures  5A-5D)  show  a  pair  of  resonances  at 

3. 2- 3. 4  MHZ  and  4.4-4. 6  MHZ  with  an  intervening  null.  The  sensitivity  to 
azimuth  is  quite  apparent.  Figure  8A  has  a  pair  of  resonances  at  3.4  and  4.8 
MHZ  which  are  probably  related  to  those  in  Figures  5A-5D.  Figure  8B  shows  a 
measured  spectrum  in  which  the  4.8  MHZ  resonances  dominates,  and  coispares 
favorably  with  Figure  5A.  Again,  the  data  show  no  pronounced  resonances  above 
6  MHZ,  for  the  same  reasons  given  previously. 

o  EALF-DOT 

The  thln-wire  results  (see  Figures  6A-6D)  show  a  single  resonance  at 

4.2- 4. 4  MHZ  and,  in  three  cases,  a  null  centered  at  2.5  MHZ  with  a  marked 
increase  up  to  the  1  MHZ  value.  This  compares  favorably  with  the  measured 
data  of  Figures  9A  and  9B.  Figure  6C,  in  which  there  is  no  direct  E-fleld 


17 


excitation  of  the  fuselage,  is  in  good  qualitative  agreement  with  the  data  of 
Figure  9C.  Note  that  the  resonance  in  the  data  occurs  at  3.8  MHZ.  The 
difference  between  this  value  and  the  4. 2-4. 4  MHZ  of  the  aodel  is  due  in  part 
to  the  absence  of  eopennage  detail  in  the  aodel. 
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Spectrum  of  JNT-DOT  for  Incidence  along  nose  obtained  from  thin-wire 


Figure  4B 


Figure  AC  Spectrum  of  JNT-DOT  for  Incidence  along  left  wing  tip  obtained  from 
thln-wlre  model. 


Spectrum  of  JNT-DOT  r  r  incidence  along  tail  obtained  from  thin-wire 


24 


Figure  5B 


Spectrum  of  JLUW-DOT  for  incidence  along  tall  obtained  from 
thin-wlre  model. 


Figure  6B  Spectrum  of  EALF-DOT  for  Incidence  et  45*  to  fuselage  obtained  from 
th  in -wire  model. 


Spectrum  of  EALF-FOT  for  incidence  along  tail  obtained  from 


Figure  7C  Measured  data  for  JNT-?VT  spectrum 


J-OOT  LUW 


238C6 . 89 


Figure  8A  Measured  data  for  JLUW-DOT  spectrum. 


I 

i 


j 


Figure  9A  Measured  data  for  EALF-DOT  spectrum. 
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D-DGT  HLF 


238C4  82 


Figure  9B  Measured  data  for  EALF-DOT  spectrum 


O-OOT  fiLF 


238C4  83 


SECTION  VI 

RECOMMENDATIONS  FOR  DATA  REDUCTION 

The  ultimate  objective  of  this  program  Is  to  provide  to  AFWAL/FIESL 
quantitative  procedures  for  removing  airframe  effects  from  the  data.  Since 
the  sensors  are  sufficiently  broad-band,  they  only  add  a  multiplicative  factor 
to  the  overall  transfer  function.  As  described  in  the  text,  the  airframe 
effects  are  analytically  described  by  the  following: 

a.  A  thin  wire  model  which  calculates  airframe  resonances,  together 
with  effective  circumference  to  obtain  surface  fields. 

b.  Frequency-independent  cylinder  and  strip  models  which  calculate 
azimuthal  currents  and  azlmuthally-dependent  fields. 

Although  airframe  resonances  are  apparent  in  the  data,  their  amplitude  is 
often  quite  small  compared  to  the  dominant  low  frequency  response.  Figure  10 
is  a  typical  integrated  waveform  for  the  surface  current  density  (not  time 
derivative)  at  JLUW.  For  such  a  waveform,  the  peak  incident  field  amplitude 
may  be  well  approximated  by  drawing  a  line  through  the  center  of  the  resonance 
and  using  the  equations  of  Table  II.  The  rate-of-rise,  however,  is  more 
strongly  influenced  by  the  resonance.  In  order  to  remove  resonance  effects 
from  rate-of-rise ,  the  following  general  procedure  may  be  used: 

1.  Determine  the  frequency,  amplitude,  and  0  of  the  resonance  from  the 
Fourier  transform.  And  check  against  the  original  waveform.  This 
requires  having  absolute  values  for  the  Fourier  spectrum.  The 
resonant  waveform  is  a  damped  sinusoid  with  undetermined  phase. 

2.  Determine  by  inspection  of  the  waveform  the  appropriate  relative 
phase  of  the  resonant  waveform.  In  Figure  10,  for  example,  the 
resonance  is  evidently  in  phase  with  the  low  frequency  response. 

3.  Subtract  the  resonant  waveform  from  the  measured  data  and  calculate 
the  resulting  rate-of-rise. 


This  method  requires  some  judgement,  particularly  in  the  choice  of 
relative  phase.  However,  for  waveforms  dominated  by  a  single  resonance,  it  is 
a  good  first-order  correction  to  the  rate-of-rise  data. 

The  brief  discussion  of  the  resonant  structure  of  the  measured  data 
points  out  the  following  information  implicit  in  the  relative  amplitudes  of 
the  airframe  resonances. 

1.  Orientation  of  the  aircraft  relative  to  the  incident  plane  wave. 

2.  The  high-frequency  roll-off  of  the  lightning  spectrum. 

Along  this  line,  a  recommendation  for  spectral  analysis  is  to  limit  the 
time  interval  to  the  leading  edge  of  the  waveform,  and  to  exclude  any 
precursor  pulses.  This  should  remove  some  of  the  "hash"  from  the  Fourier 
transforms . 
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